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ABSTRACT 


A review of the current theory and methods for the computa- 
tion of maximum flow in networks is presented along with a simplified 
method for determination of a feasible flow in networks with upper 
and lower bounded arcs. A computational procedure is presented which 
is used to calculate the maximum flow for a general network. The 
network is reduced to an equivalent basic network. An associated 
network is used to compute a feasible, then the maximum flow for 
the basic network. A computer program is included for use in 


computation of maximal flows in large networks. 


DUDLEY KNOX LIBRARY 


— NAVAL POSTGRADUATE SCHOOL 
DUDLEY KOO ADUATE SCHOOL MONTEREY CA 93943-5101 


MONTEREY CA 93943-5101 
TABLE OF CONTENTS 


Chapter Title Page 
I. Introduction 7 
II. Review of the Literature 1] 


III. A Simplified Method for Computation of Maximal Flow 17 


IV. Procedure for the Solution to the Generalized 


Maximum Flow Problem 20 
Appendices 

A. Detailed Flow Charts 34 

B. Fortran Listing 48 





LIST OF ILLUSTRATIONS 


Figure Title Page 
Ae: Generalized Network G 24 


2. Generalized Network With Single Arcs Replacing 


Multiple Arcs 24 
o. Equivalent Basic Network 26 
4. Transformed Network G' 26 
General Flowchart 30 
Detailed Flowchart 3/7 





CHAPTER I 


INTRODUCTION 

The growing uses for mathematical programming models in trans- 
portation and communication networks have prompted study in the devel- 
opment of computational methods to determine maximal network flows. 
The subject of this thesis is the mathematical problem of determining 
maximal steady state flows in networks which are subject to capacity 
limitations on the arcs and nodes of the network. 

An efficient and widely used method to compute maximal flow 
was developed by Ford and Fulkerson [7]. This method along with a 
review of the current theory will be presented in Chapter II in 
order to give the reader the theoretical foundation upon which the 
remainder of the thesis is based. Chapter III introduces a new 
method which utilizes the Ford and Fulkerson algorithm to compute the 
maximal flow in networks with non-zero lower bound limitations on the 
flow in the arcs. Since many networks of interest are more general 
in nature, the generalized network is introduced in Chapter IV along 
with a computational procedure for the determination of maximal flow. 
To illustrate the techniques proposed, an example is included. The 
procedure can be used to solve small problems by hand, but the amount 
of work increases rapidly with the size of the network. Therefore, a 
computer program which may be used to solve larger problems is 
included in the appendices. A brief description of the program is 


found in Chapter IV. 


Definitions and Symbols 


In order to establish a common understanding as to the exact 
meaning of various terms and symbols used to present the material of 
this thesis, an initial set of definitions is presented. Other terms 
will be defined throughout the material as needed to facilitate the 
presentation. 

Define a network, G, as a graph which can be represented on 
a plane in such a way that the set, N, of nodes are distinct points 
and the set, A, of arcs are simple curves which connect two distinct 
nodes. Furthermore, no two arcs can meet except at the nodes which 
are their extremities. The graph contains no loops. 

Define w (N) to be the set of arcs or flows incident to a set 
of nodes, N, and w(N) to be the set of arcs or flows incident from N. 

A directed arc, aay is defined as an arc incident from a node 
i and incident to a node j. 

An undirected arc connects node i and node j without having an 
orientation. 

Flow can be thought of as a value of the steady state rate of 
movement of a homogeneous commodity along a path or channel. Positive 
flow in a directed arc, ayy? will move with the orientation of the 
arc and will be denoted by the symbol, Yay: Positive flow in an 
undirected arc, denoted Yay may move in either direction but not 
in opposite directions at the same time. 

To each arc, ay c A, there will be associated two real 


numbers, b. ., and C, 5? which represent, respectfully, the minimum 


Lb 
and maximum allowable flow in that arc. The bounds associated with 


undirected arcs, ai 4? are denoted Dy oy? and Ci 4: 


An arc is saturated if Ys = Cig: 

Define Y as the set of flows in a network G. 

Initially, the capacities of the nodes will be assumed to be 
infinite. The nodes will be classified according to the value of the 
difference of the sum of the flows into the node minus the value of 
the sum of the flows out of the node. A node is designated as a 
source node if this difference is negative, and as a sink node if this 
difference is positive. For all intermediary nodes, the sum of the 
flows into a node must equal the sum of the flows out. The set of 
source, intermediary and sink nodes will be denoted as S, I, and T 
respectfully. 

Define F as the value of the flow in a network. 


F = a 
J; it (T 
A flow is feasible if and only if: 


A 
1,t 


aes 27 ),5 SP ye 


(2) F- y y = 1) 
Vg, gow (8) § 
> aE B40 
Y4 4S (ty tot 
(3) -0 4§ 61 


y, you" (1) Se yy yout(g) 1 
Define a chain as a sequence of arcs such that each arc, ayy 

i,j ¢ I, is connected to an arc ay 4 at node i and to arc ay at node j. 
A cycle is a sequence of arcs where: 


(1) The sequence is a chain. 
(2) The sequence does not use the same arc twice. 
(3) The initial and terminal nodes of the chain coincide. 


A cut is defined as any set of directed arcs containing at 
least one arc from every chain of positive capacity joining the 
source to the sink. The value of the cut is equal to the sum of the 


capacities of the arcs of the cuts. 


Basic Network 

Much of the initial development of the material in this thesis 
deals with computation of flows in simple or basic networks. There- 
fore, for the purpose of this paper, a basic network, G, 1s defined 
as a network which has the following properties. 

1. The network has a return arc, and? which is added to the 
network for computational purposes only. 

2. Node 7 has only one arc, an? incident to it. Node 1 fs 
designated as the source. In the basic network there jis 
only one source node. 

3. Node n has only one arc, an? incident from it. Node n 
is designated as the sink. In the basic network there is 
only one sink node. 

4, For each arc, ayy the associated Cy 4 is positive and 
by = 0, 

5. Two distinct nodes (i,j) may be connected by only one 


directed arc. The network contains no undirected arcs. 


CHAPTER II 


REVIEW OF THE LITERATURE 
In order to acquaint the reader with the subject of maximal 
flow in networks, a review of the current theory and primary combina- 


tional techniques is presented in this chapter. 


Basic Network Theory 

A great deal of theoretical work in networks has already been 
done. A few of the relevant theorems are presented here in order to 
provide justification and a rationale for the combinatorial techniques 
and algorithms that are presented in later sections. 

Theorem 1. Opposite directed flows on the same arc cancel [4]. 
This theorem means that given an undirected arc, ayy? with flows 
Yi > 0 and Ys > 0, these two flows can be replaced by the flows 


| 
‘abba lesic: 


yy 

Theorem 2. A set Y with F>0, satisfying the capacity con- 
Straints and node conservation equations can be decomposed into a set 
of positive chain flows from the source to the sink and a set of 
circular flows such that the direction of positive flows in any 
common arc is the same for all chains [4]. 

Theorem 3. There exists a positive flow from the source to 
the sink if there exists at least one chain of arcs with positive 
capacity from the source to the sink [4]. 

Theorem 4. A flow F° is maximal if and only if the maximal 
flow is zero in a second network formed by replacing Ci Dy c.f 


TJ 


ee 
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Theorem 5. Given any partition of nodes into two sets, where 
the first set includes the source node and the second set the sink, 
then a feasible solution F is maximal, if every arc, aay that joins 
a node in the first set to a node of the second is saturated [4]. 

Theorem 6. The maximum flow in a network is equal to the 


minimum cut. 


Tree Method 

Dantzig [4] has developed a systematic procedure for com- 
puting the maximal flow in a network by using chains connecting the 
Source and the sink of the network. A knowledge of the tree method 
will give a better understanding of the Ford and Fulkerson algorithm 
which is presented in the next section because the two methods are 
identical in principle. 

A tree [1] is defined as a connected graph with at least two 
nodes and no cycles. 

Consider a basic network such that: 


iden eS 


> 0 Var 7c 

The arcs may be directed or undirected. Initially all arcs 
are unmarked and the feasible flow is zero. In order to determine 
the maximum flow: 
1. Create a tree such that the arcs are unsaturated. 
2. Select two sub trees - one T,> branching out from the source, 
node 1, and the other, Ty» branching out from the sink, n, such that 
each node is reached by just one arc of the tree. 
3. Choose one unsaturated arc which connects the two trees. Thus, 


there will be just one chain from node 1] to node n. If none can be 


found, go to step 7. | 


TZ 


4. Determine the amount that the flow can be increased along the 
chain from node 1 to node n, by the following procedure. Determine 
the amount that the flow can be increased, AX, where AX = min (Cc; - 
Yay Ye y): Flow along the chain directed toward the sink is 
designated Yi and flow directed toward the source along the chain 
is designated Yer There are two cases: 
(a) AX > 0, increase the flow by AX along the chain, go to 
next step. 
(b) AX = 0. This means one arc in the chain is saturated 
by a flow directed toward the sink, go to next step. 
5. Mark one of the saturated arcs. 
6. Eliminate the saturated arc from the chain. Thus, there are two 
sub-trees. Return to step 3. 
7. Check all marked arcs joining Ty to Th If all these arcs have 
flow from Ty to i then the optimum has been reached. If any arc 
has flow from Te to T)> use this arc to connect Ty to Ti, and proceed 


with step 4. 


Ford and Fulkerson Algorithm 


Given the basic network, the most direct method developed to 
date of computing the maximum flow is the algorithm developed by Ford 
and Fulkerson. This method is the one most cited in the literature 
concerning maximal network flow. The procedure is outlined below. 

1. Begin with a feasible flow in 0 in all arcs. 

2. Label each node as follows: 


a. Label node 1, with the number 1. 
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b. If node i is labeled and node j is not labeled, label 
node j with the number (i) if: 


(1) as exists and Yi 5 < SHG 


(2) as exists and “sg > 0 
c. If the sink, node n, becomes labeled by this procedure, 
then the flow from node 1 to node n can be increased. 
3. If node n is labeled, construct a simple chain from node | 
to node n by back tracking from n to 1 using the labels 
on the nodes. 


a. The chain created will have the property that all arcs, 


a. j on the chain directed toward the sink have a flow 


O<y. and all arcs a, , directed toward the 


nace 

a TsJ 
source will have a flow Ya? Oy: 

b. The flow can be increased by an amount AX where: 
AX = min(c 


~ y 3 ome Y sees ay yk Nn chain. 


15J 15J 
4. Increase the flow along the chain by an amount AX and 
return to step 2. 
5. This procedure is repeated until the sink cannot be 
labeled. Flow is maximized when this is the case. 
Now consider the more general problem in which the lower bounds 
on the directed arcs are no longer assumed to be zero. Thus the flow 
in each arc is bounded both above and below. 


Such that: - <b, < + a A 


15-715 = S15 = ij 
Theorem 7: In the case of networks with lower bounds on the 
arcs, C is defined as the subset of N such that if node 1 C the node 


n 4 C or if node 1 c C then node nc C. 
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A flow y; j exists and is feasible if and only 1f: 


r. 


Pigg ts, Ogre” 
Cc > x b 
a Ly ml + T5J 
a; 3 fw (C) aa cw (C) 


This more general type of problem can be satisfactorily worked 


by using the Ford and Fulkerson algorithm if a feasible flow can be 


found. 


The Ford and Fulkerson algorithm can be used to solve for a 


maximum flow in an associated network G' which is derived from the 


original network, G. This method is outlined below. 


1. 


Create an associated network G' such that G' contains 
the nodes and arcs of the original network G with arc 
capacities C43 = Ci : bs . The lower bound, 
big equals zero. 
Add a new source, labeled 0, and a new sink, labeled m, 
ijn + 1), to G’. 

a. If the arc, aay? exists, then construct two new 


arcs according to the value of the lower bound, b. 


oJ 
If b. j > 0, then construct the arcs, a’ a of capacity 
' = ‘ : : Q = 
Cos ym bs and a 0,4 with capacity c 0,4 by : 


If b. . < OQ, then construct the arcs, a'. ; With 
15J 0,1 


Capacity c . b. 


r j and a 5, with capacity 


ee 
b. The Ford and Fulkerson algorithm can then be used to 
determine the maximum flow in G' from node 0 to node m. 
Two cases can result from step 2b. 


a. Case one: If the flow in each arc, a 0,3 and a 5m 


ifs 


j c N', equals the capacity of each arc, then the flow 
is feasible for the network G and the problem has a 
solution. 

b. Case two: If for any arc, a 0,4 or a 34m 
the flow in that arc does not equal the capacity of 


jcAt, 


the arc, the problem has no solution. 
If a feasible flow for G can be found by the above pro- 
cedure, transfer the flow from G' to G by using the 
following transformation: 
Yig = Jae i bs V aa cA 
Ford and Fulkerson's algorithm is then used to determine 


the maximum flow using the feasible flow as a starting 


basis. 
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CHAPTER III 


A SIMPLIFIED METHOD FOR COMPUTATION OF MAXIMAL FLOW 
To compute a maximum flow in networks with non-zero lower 
bounded arcs, an initial feasible flow must be determined. The method 
presented in the last chapter produces a complex associated network, 
G', even with a relatively simple network, G. The principle advantage 
of the method presented here [10] is that an initial feasible flow 
can be determined by use of a more simple network, G', than in the 


previous method. 


Proposed Algorithm 


The development of the algorithm is based upon a consequence 
of the flow balance equations. In order to provide a rationale for 
the procedures presented, the development of the associated network 
will be presented now. 

Given a basic network G where the flow in each arc is bounded 
as follows: 


co -< f) 


ee an ee oe oe Sa ager Cah 
Sj ai i,j- 


The problem of determining the maximum flow in G can be stated as 


follows: 
Maximize: | 
Subject to: C Yea . y¥. .=0 View 
a Nonds | ty see lye 
aaa © (i) Yi jet (i) 
Bigg SOS Gre rag 


Now transform the entire network G into an associated network 
G' using the change of variables: 


A. Sa ee eye ee ». 
Yi j Y 4,3 Daj a oe 


the problem then becomes: 


Maximize: y ae 
Subject to: 
\ ae y' b 
' “fe Jol t +). T»5J C dual 
y',.3 ge MW) y';3 ce MW) bs w (7) 
= YicN 
+). Ts5J a 
by Curt i) 
ae eet, Cageemeieg 
Let . a i cN 
-72\ joi + ree | 1 = 
b, 7 cw (i) b,j cu (i) 


Note two things. First the sum, 2 d., is equal to zero. 
Second, if a flow of -d. was passed ehrolifit each node, the problem 
would be exactly the problem of flow maximization in a basic network. 
Hence, a technique similar to Berge's [2] can be used in which new 
nodes 0 and m = ntl are added to the network G'. Therefore, only a 


maximum of one arc per node is added to G' in this method. 


Procedure 

Using the results of the previous section the solution to the 
maximal flow in the above network can be determined using the following 
procedure. 

1. Add a new source and sink node, denoted 0 and m (m=n+ 1), 
sampect Ful ly TOMmGns 

2. The arcs a8 ‘ and a; m are added as follows. If d. =. 0), 


thevan arc, ag a> is constructed with capacity d.. If d. < 0, an 


arc, a" m? is constructed with capacity -d.. If d. = Q, no arc is 


added to G' for node i. 
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3. The Ford and Fulkerson algorithm is then used to maximize 
flow in G' from node 0 to node m. Again as in the previous method 
this maximum flow is equivalent to a feasible flow in the network G 
if and only if: 

=O. Yiad>Q0O 
¥Y13ada<0 


te 
i} 
Co 


4. If the problem is solvable continue by eliminating all 
arcs incident from node 0 and incident to node m. Node 1 and node n 
are now the source and sink respectfully for G'. 

5. Using the flow generated in step 3 as a starting basis, 
maximize flow from node n to node 1 by again using the Ford and 
Fulkerson Algorithm. 

6. Upon completion of step 5, transfer flow from G' to G by 
use of the transformation equation: 


y ec is A Scone Yiaee. clA 


La T3J Us 


The maximum flow in G is equal to vm 


CHAPTER IV 


PROCEDURE FOR THE SOLUTION TO THE GENERALIZED MAXIMUM FLOW PROBLEM 

A proposed methology for the computation of a solution to the 
generalized maximal flow type problem as defined below is presented. 
The generalized network is first reduced to a basic network and then 
the technique presented in Chapter III is used to solve for the equiv- 
Nalent maximum flow in the basic network. An example is presented in 
order to demonstrate the techniques discussed in this chapter. A 
description of the computer program contained in Appendix A for use 


in solving maximal flow problems for large networks is included. 


The Generalized Network 
Generally maximal flow solutions are required for networks 
which do not fit the restrictive definition of a basic network. For 
the purpose of this thesis the generalized network is defined as a 
modified basic network with a set of arcs A, and a set of nodes N, 
such that: 
1. Each set w (i) and w' (i) can contain any number of 
directed and undirected arcs. 
2. Let D be a subset of A such that if and only if the arc 
ys is directed it is a member of D. 
3. Let U be a subset of A such that if and only if the arc 
as is undirected it is a member of U. 
4. K is a subset of N such that if the node i is bounded 


it is a member of K. To the bounded node, i, we associate 


two real numbers, bss the lower bound, and, C., the upper 
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bound such that the total flow through the node i, 
denoted Yas is bounded: 
b. Sh tse. 
5. All nodes, j, not bounded belong to the subset 1, where 
ec N. 
6. The network may or may not contain multiple sources or 


sinks, but shal! contain at least one source and one sink. 


Procedure for Reduction of the Generalized Flow Problem 

In order to solve for the maximal flow in the generalized net- 
work by the techniques presented in Chapter III, the network must be 
transformed into an equivalent basic network. 

The proposed methods for reduction are based upon results of 
the general theorems presented in Chapter II. References cited give 
further proof of some of the methods used. 

First, consider the case in which the generalized network, G, 
contains both directed and undirected arcs. This network can be 
simplified by replacing each undirected arc by two directed arcs as 
follows: 

Given an undirected arc a. ; cl 


3 


Replace a. 


‘sj by agg and a... c D 


TJ JI 
Where 


A 


ee — es 


= jqimiets,i ~ Sing 
Upon reaching a solution for the maximal flow, the orientation 
by the flow in the undirected arcs can immediately be determined from 


the values of Ys and y; j Since one flow will be zero. 


j 


2] 


Now consider a network in which two distinct nodes, i, j, are 
connected by more than one arc. Successive arcs, flows, and bounds 


are denoted by prime symbols. For example if there are m arcs from 


th 


node i to node j, the k~ arc would be denoted a. a It will be 


assumed that each arc has both upper and lower bounds. 
(k) (k) (k) = 
bs Ses a k = 1,2,...m 
The set of arcs a. Riad k = (1,...m) may be replaced by a single arc 


: z (k) Da, oh) 
ij) oe 1.  ) 


Next assume a pair of nodes (i,j) are connected by two oppositely 
directed arcs aa and as ee 
The pair of arcs may be replaced by a single directed arc 


in either direction as follows: Ifc. .>b.. 


TsJ ~ Jol 
Replace a and as by aa such that:° 
! = = 
ae 
Do aD - C 


isj Ted lid 

Otherwise orient the arc from j to i. 

If a network has multiple sources and, or sinks, merely create 
a single source and, or single sink. Then construct one directed arc 
of infinite capacity from the new source to each old source and one 
directed arc of infinite capacity from each sink to the new sink [8]. 

If the flow through the nodes is restricted to lie between 
two values, the character of the maximum flow problem itself remains 
unaltered. Assume we are given a network in which the flow is subject 
to both arc and node capacities, where all arcs are directed or have 
been converted to directed arcs by the previous techniques. This 
network can be reduced to the basic network maximum flow problem by 


the following procedure. Expand the network such that for each bounded 
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1\ 


node (i) c K, two new unbounded nodes (i', #'') r2 created, so that 


LL joins the two new nodes with the 


upper and lower boun.is associated with the bounded node in the reduced 


the single directed arc, a 


network. To each arc, ap as create an arc a’ Ar in G' and to each 
arc a. : in G create an arc Aa ; in G'. The bounds on the other 
b ] 9 


arcs remain the same. 


Example Using the Simplified Technique 


Given the g@ucralized network show? in Figure 1, we are 
required to compute the maximal flow from node 2 and 3 to node 10. Or; 


Maximize F = y Bye = y oy Ty, 
Vine (ny in Sad 7240 6,10 


The flow in each arc is bounded as follows: 


aer< oo 5 ] <%657 ES ocd) 3 < %6 105 9 

] 2) 
en, 2 * PEWS <> 2 
4 ane |¢ 


The flow in all nodes except 5 and 6 is unrestricted. The flow in 
node 5,6 is restricted as follows, 


1 eels 4 < ye < 10 


This flow problem can be simplified by using the techniques 
presented in the last section. The procedure is outlined by ordered 


steps. Although many steps could be combined for this simple 
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FIGURE 1 


GENERALIZED NETWORK G WITH UPPER 
AND LOWER BOUNDS ON ARCS AND NODES 
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PIGURE 2 


GENERALIZED NETWORK WITH SINGLE 
ARCS REPLACING MULTIPLE ARCS 
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example, it is felt that for hand and computer computation this is the 
most efficient procedure. 
1. Add the single source node, denoted as 1, to the network. 
In addition, add the arcs ay 2 and a1 3 to the network 
such that: 
O< V9 Sale 
0<¥,3<+* 


2. Replace each undirected arc a. j with upper bound c, . > 0 


15J 
by two directed arcs a. . and a. 


15J Js1° 
Therefore: 
- (2) (2) (2) (2) 
an 5 becomes ay 5 and ara 0 Lae < 5598 Sena eee 
: (3) (3) 
ay 6 becomes an 6 and a6 4 0 < YA 6 Kaci) 7 ¥ SMC 


a5 6 becomes a5 6 and a6 5 0 5 Ye 6 ws, 0 =I 6s <0 


3. Combine all arcs joining each pair of arcs (i,j) according 
to the following procedure. 


a. Combine all arcs which are similarly directed, into a 


Single arc. Thus the arcs a al?) a3) become an 6 F 
00 rr re ae = 
“1,6. 6 “ho meee 7 C ? 
(URE (2) , gle | - 
be Fl wae tee | tC OPO 


Similarly form the arcs an .5° 85 4 and a7 6 23 


Cas = 13 Das = 3 
Ce 4 = 13 bs 4 = 3 
C7 ¢ = 8 b7 ¢ = 8 
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FIGURE 3 
EQUIVALENT BASIC NETWORK 





FIGURE 4 


TRANSFORMED NETWORK G' 
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b. Transform each pair of arcs which are directed in 
opposite directions between each pair of nodes 


into a single arc. 


Thus: an 5 and a5 4 become aq 5 23 
Cas = Cas - bs 4 = 13-3 = 10 
ba 5 = ba 5 - Co 4 = 3-13 =- 10 


Similarly the pair of arcs an 6 and a6 4? a5 6 and 
A653 85 7 and a7 53 87 ¢ and ag 7 are combined. The 
simplified network is shown in Figure 2. 

Complete the transformation to the basic network by the 

following procedure. 

a. Expand the network so that all nodes have infinite 
capacity for flow as outlined earlier. Thus for 
nodes 5, and 6 the new arcs a5 9g and dg g are 
created and the equivalent network shown in Figure 3 
is established. 

b. Add the return arc 410.1 for computational purposes. 


Transform the network G to G' by using the change of 


variables equation yee ve ‘i Dy aa cA 
Hoy = Y 10M pea We v1) .3 72 
on - Y oe &) Yo.6 we 2.6 Vaui gomes 
Vag Soaygeete Ya,5 8% a5 - Wedge 8 ayer | 
leper X50" Yes = ¥'c Jo + °F%7,6 7 Y 746 
ae aio’ ° Yee “Yoo,c - "See Yoel 


2 


vemion 3.10 ° © Yoo ~ Y 9:10 


2] 


The equivalent problem is 


obtained: 


Maximize Y'10 1 with the node flow satisfying the balance 


equations as follows: 


Node 
Pon) eae, 1,2 = 0 
ME 3) OREN AE = 
oe) 34 ees go 
a! eee De ie emea.6 +o “4 
SY gage 9) a5 Gee * “a 
em ee em eee Ge 
1 Vege) 7.6 9) 7a a 

Jone. 116.6  Saowmeer io | a 
a ee noe | = 
i) eai0 Meio ace memo 3 
Therefore the values of the remainder d. are immediately 
seen to be as follows: 

d,=0 d,=-8 d=-8 dy = + 22 

d,. = - 9 de = - 13 d_ = - 5 de = +1] do = +] dig =+9 


Using the values of d., i 


new source, sink and arcs 


= 1,,,,10, as the capacities, the 


are added to obtain the equivalent 


network G' shown in Figure 4. 


Then using the Ford Fulkerson algorithm to obtain a feasible 


flow in G the following values are obtained for the flow in 


Ge 
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ieee eee Yo * 8 YG Get 
y3,4 © 9 yes ee Wie? Yas = 8 
y58 1 / oem” 2G, ~ ° Yano * 
13,6 2. mwa! Yeo” ° YS,105 ° 
y10,0 ee 


6. Note that the flow in all] sink arcs and source arcs is at 
arc capacity, therefore, a feasible flow exists. Dropping 
nodes 0 and 11 plus the associated arcs, maximize flow 
from node 1 to node 10 using the Ford Fulkerson algorithm. 
In this case the feasible flow computed is the maximum 
flow. Then, using the change of variables equations the 
resulting maximal flow on G is determined. Finally G 
is simplified by eliminating arcs a5 g and 46 9° The 


maximal flow is 20. The flow in each arc is as follows: 


4 | Oa et °° abe 
11,5 | peteee’s.4 = > 5,7 3? 
Yougig “ene. 5  ° eg 

6 


Ye .10 = Mame 5 = > 970 = 


The flow in node 5 is 9 and in node 6 is 9. 


Computer Program Description 


Since it is often the case that it is necessary to solve the 
maximal flow in a much larger and more complicated network, a computer 
program for the solution of the generalized maximal flow problem as 
presented in this thesis has been written. The program is written in 


FORTRAN IV for use on an IBM 360 computer. 
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START 
READ IN 
G 


ALL ARCS INCIDENT 
TO AND INCIDENT 
FROM THE SAME SET 


OF NODES ARE 
COMBINED 





NUMERICAL VALUE 

OF THE TERMINAL 
NODE IS INCREASED 
TO PERMIT EXPAN- 
SION OF THE NETWORK 


COMBINE MULTIPLE 

ARCS JOINING TWO 

NODES IN OPPOSITE 
DIRECTIONS 


TRANSLATE BOUNDED 
NODES INTO 
BOUNDED ARCS 


COMPUTE CAPACITIE$ 


SINK NODES IN G' 





CREATE G' 


GENERAL FLOWCHART 
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f 
i 
t 
| 


ere 9 ee, Brees Ae ee ee eee 


| 
| 
: 







CREATE CHAIN 
FROM SOURCE TO 


SINK 
LOAD PATH . Bae 
dene << LABELED 
es . 
COMPUTE MAXIMUM THR 
FLOW INCREASE FOR 
| G'? 
ae ae YES 
FLOW IS INCREASED! 
ALONG PATH 
| 
i ELIMINATE ALL 
ARCS FROM SOURCE L. YES FLOW 
AND SINK NODES FEASIBLE 
OF G' ? 
NO 
CLEAR 
NO FINAL 
FLOW NETWORK 
ARRAY ? 


GENERAL FLOWCHART CONTINUED 


3) 










TRANSFORM FLOW 


FROM G' TO G 





TRANSFORM BOUNDED 
ARCS TO BOUNDED 
NODES 












SEPARATE APPRO-; 
PRIATE ARCS AN 
FLOWS TO TWO 
OPPOSITE ARCS 
AND FLOWS 













er 5 


CHANGE VALUE OF 

TERMINAL NODE AN 
COMPUTE MAXIMUM 
FLOW FOR G' 


PRINT ARC FLOWS 
AND MAXIMUM 
FLOW FOR G 


The maximal flow in generalized networks is computed using 
the methodology outlined in this chapter. In addition, a set of 
routines is added to the program which will permit the user to 
transfer several networks to input data directly without first 
reducing the networks to basic networks by hand or auxiliary methods. 
However, as the program is written each network must have a single 
source and sink and all undirected arcs are to be represented by 
two oppositely directed arcs. The program has been checked out and 
functions properly. 

The present program is written for a 50 node network. In 
order to solve larger networks the size of the arrays in the dimension 
Statement must be increased. The general flow chart of the entire 
program is given on pages 30 and 31. For the interested reader a 
detailed flowchart is presented in Appendix A. A Fortran listing 
of the program is contained in Appendix B along with the solution 


to the previous example. 


Conclusions and Recommendations 

The general maximal flow problem can be solved more efficiently 
by use of the method developed by H. Greenberg. ‘The advantages of the 
technique are evident when solving small problems by hand. The tech- 
niques presented provide a systematic approach to obtaining a solution. 
The problem of computing the maximal flow in a mixed network in which 
the flow in the undirected arcs are bounded below by a positive 
number cannot be solved by this method. In fact, to the best of my 
knowledge this problem is unsolved in the general case except by 


complete enumeration techniques. 
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APPENDIX A 


DETAILED FLOW CHARTS 
This appendix contains the detailed flow chart of the computer 
program in addition to a description of variables used. The statements 
are written in FORTRAN IV with arrows indicating the logical flow 


through the system. The input and output statements are not enclosed. 


Description of Variables 


Variable Name Description 

Gil sa) ) Capacity of the directed arc from node I to node J. 

Bt eel) Lower bound of the directed arc from I to node J. 

Y(I,J) Initially used to read in data. 

Y(I,1) The number of the node the arc as is incident 
from. oJ 

Yells. 2.) The number of the node the arc a. : is incident to. 

VGiss) Upper bound of arc a. 

Y(1,4) Lower bound of arc a. F 

Y(I,J) The value of the flow in a directed arc from I 
to J after input data is transmitted to arrays 
C(I,J) and B(I,J). 

CEC) Dummy array used to compute values of the 
capacities of the source and sink arcs for G' 
and then the label of the nodes in the label 
routine. 

NL(I) Dummy array used in the label routine as a check 


to determine if all nodes have been checked and 
then to store the path used to increase flow. 


XNC(1,1) Number of the node that is bounded. 
XNC(1,2) Upper bound on a flow through a node. 
XNC(1,3) Lower bound on flow through a node. 
XNC(1,4) Value of flow through a node. 


NN Number of networks. 
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N Number of nodes in the original problem. 


M Value of the designated number of the sink node. 
KX Value of the designated number of the source node. 
NA Number of arcs in the original problem. 

NC Number of bounded nodes. 


Input Card Format 


All entries on the card format in this section are right, 
justified and are fixed point entries if the variable name begins 
with I, J, K, L, M, or N and are floating point entries for all other 


variable names. 


Card #1 
Column Name 
1-4 NN 
Card #2 
1-4 N 
4 - 8 NA 
8 - 12 NC 


Card #3 through Card NA + 1 


1) ano Y(1,1) 
I= 20 ¥GeZ) 
2) =ee0 Y¥(1,3) 
31 a0 ¥(1,4) 
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Card #NA + 3 through Card NC + 3 





Column Name 

1 - 10 XMC(I,1) 
11 - 20 XNC(I,2) 
21 - 30 XNC(I,3) 


Detailed Flowchart 
The detailed flowchart of the program is given in the next 


1] pages. 
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=e 


B( id) = 0 
C(I,J) = 0 


DO tf = 1,NA 


Gu tte me 


oe 


= ¥(T,1) + 1 
= ¥(I,2) + 1 





DETAILED FLOWCHART 


uy 


L=N+1] 
K=N+NC + 1 


DO I = 1,M 


Cink) = Colas 
B(I,K) = B(I,L) 


| 
| 
: 
| 





C(J,1) = - C(J,1) 


C(I,J) = -C(I,J) 


DETAILED FLOWCHART 
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ie 
: NO 
< NC > 0 





DO I = 1,NC 


— = XNC(I,1) + 1 


Ne” CA ital owas |_TSS 


‘Ss 7 
a 


C(L,J) = 


! 
~O 
— 
~~ 

w 
Co 
— 


B(L,J) = B(K,J) 





DETAILED FLOWCHART 
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et ay 


Pe OPE At, a FE OE = 


LT A TP Meal NT TDI LO ERRAT AN ver ae oD Lr 3 WEEE all os 


SSS er Sy SE EE Tt ee eo i | ge 





pales et ay POOL © -06-t al Mini ctype 


YES NO 

| 

| 

, 

| C(1 5d) = Cuea)) YES 
: =B(Td) 

: No 

: 

| 

} 


SO = SO + B(I,J) } 








DETAILED FLOWCHART 
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C(M-1,2) = .9999 


DO I = 2,Mx 






C(I,M) = -xL(1) |_NO VES 1 (1,1) = XLQ) 





DETAILED FLOWCHART 


4] 


XL(KX) = KX 








DO K = 1,LN 
ES 


DO J = KX,M 













. 
3 ! NO Y 
_ 
: | NO 

| 

| 

| C(J,1) > 0 

| Le 

: 

j 





NO NO 


Yes 


a eee 


DETAILED FLOWCHART 


4? 


I=] 
NL(I) = M 

K = NL(I) 

T=1+ 1] 
NL(I) = XL(K) 

K = NL(I) 


YES 


X = 9999.9 











= 1,M 
NL(I) 
NL(1I+1) 


VES 


X = C(L,K)-Y¥(L,K) 
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. 
| 
© 

T- KR es 

ll | tl 

— ao = 

rer « 

—_—o— 

— me 

— te 


ERO OUT ARRAYS NO 





Ee 0 
NL AND XL 
ra 
YES 
Y(L,K) = NO C(K,L) > 0 YES V(K,L) = 
Y(L,K) + X 
N Vike 
=a y 
oD YES =-N+NC +1 we 
=e = 
NO 
MJ =M- 13 





KX = 2 
M=M- 1 
C(M,KX) = 0 


wy 
DETAILED FLOWCHART 
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Pane 






NO 
XNC(1,4)+1 





rave ey ee eos’ 


DO J = KX,M 





* Sl ee ee ee A ee ee ee ee 
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DO J = KX,M 






Y(J,1) = Bien) 
¥(1,0) = Y(I,J) 
+ B(J,I) 







Biles +. C(alsse) 





hon = B(T,J) 
- Y(I,J) 


Y(I1,J) = B(I,J) 









L=N+1 





K = NotaNC +31 






DO I = KX,M 
C(I,L)= C(I,K) 
B(I,L)= B(I,K) 


V(I,L) = ¥(I,K) 
C(I,k) = 0 






DETAILED FLOWCHART 
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PRINT 
SOLUTION 


OF NETWORK 
FLOW 


ALL 
NETWORK 
OWS COMPLETEL 












ZERO OUT 
Y -ARRAY 








RETURN TO 
START READ 
IN NEXT 
NETWORK 







DETAILED FLOWCHART 
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APPENDIX B 


FORTRAN - IV LISTING 
This appendix contains the complete listing of the computer 
program. The problem illustrated in Chapter IV was used to display 
the output format. Comment cards are included for clarification and 


identification of the main sections of the program. 
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